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Abstract. We investigate the exact solution of the q-deformed one- dimensional Bose gas to derive all 
integrals of motion and their corresponding eigenvalues. As an application, the thermodynamics is given 
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1 Introduction 

One-dimensional Bethe Ansatz (BA) solvable lattice mod- 
els are of intrinsic theoretical interest. On the one hand, 
the BA solution allows it to diagonalize the Hamiltonian 
and - if the number of lattice sites is L - a set of further 
L — I independent operators, the higher integrals of mo- 
tion. As an application, ground-state properties as well 
as the thermodynamical behavior of the model may be 
obtained exactly. On the other hand, much insight can of- 
ten be gained by formulating an effective field theory for 
low energies - i.e. for long-wavelength excitations. If the 
low-lying excitations are gapless, such a continuum the- 
ory is relativistically invariant and displays the universal 
properties of the model. 

In this work, we investigate the one-dimensional q- 
deformed Bose model, which can be viewed as a lattice 
regularization of the one-dimensional Bose gas. The ex- 
act solution of this (7-Bose model was found as a special 
case of a more involved g-deformed quantum lattice-model 
in [4,5], and obtained directly via the algebraic BA in 
[9]. It turned out that for small coupling, the continuum 
limit of this model is the one-dimensional Bose gas with 
^-interaction, also solved by BA [18]. This model can be 
cast into the form of a Luttinger liquid [14,11]. The cor- 
responding critical behaviour of the g-Boson model was 
analyzed in [6]. Very recently [12], the Luttinger liquid 
fixed point at zero temperature together with the integra- 
bility of the model were used to calculate the three-body 
correlation function. 

Our motivation to reconsider the g-Boson model is 
twofold. First of all, we would like to point out the di- 
agonalization of all integrals of motion (IM) in the frame- 
work of the algebraic BA. In this approach, the IM are 
obtained from a generating functional, the transfer ma- 
trix. We use a slightly different, but equivalent, trans- 
fer matrix compared to the one employed by the authors 



cited above. This allows us to obtain rather simple ex- 
pressions for the IM. Our second motivation is the study 
of thermodynamic quantities like the specific heat C(T) 
and the charge susceptibility Xc{T)- Although the equa- 
tion for the thermodynamical potential has already been 
found in [9] , no attention has been paid yet to the actual 
calculation of (experimentally accessible) susceptibilities. 
Apart from the physical insight that these quantities of- 
fer, their low-temperature behaviour is of special interest: 
it reveals the operator-content of the effective field the- 
ory including irrelevant operators which necessarily ap- 
pear when performing the continuum limit. Whereas it 
is due to the Luttinger fixed point that limT^o C(T)/T 
is constant, the leading T-dependent corrections to these 
constants are governed by the leading irrelevant opera- 
tors. The same is true for a low-temperature expansion 
of XciT), as specified in section 3. For the XXZ Heisen- 
berg chain, these were determined in [19] and excellent 
agreement with the finite-temperature magnetic suscepti- 
bility of that model was found in [21]. In this work, the 
leading irrelevant operators of the q-Bose model are given, 
without determining their amplitudes. 

This article is organized as follows. In the next section, 
we introduce the model, summarize its exact solution and 
give the IM. The third section consists of the formulation 
of the effective field theory, whereas the fourth section is 
dedicated to calculating both ground-state and thermo- 
dynamic properties from the exact solution, while making 
contact with the field-theoretical predictions. We end with 
a short summary. 



2 Integrability of the g-Bose model 

In this section, we introduce the model and summarize its 
exact solution. 
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2.1 Algebra of observables, Hamiltonian and Baxter 
equation. 

Let b, b^, n be the generators^ of g-oscillator algebra: 

bbt = 1 - g2n+2 ^ btb = 1 - , 

bn = (n + l)b , b^n = (n - l)h^ 

with < g < 1. The g-oscillator is represented in the Fock 
space J^: 

n|n)=n|n); n = 0, 1, 2, 3, . . . ; b|0) = . (1) 

Consider a chain of length L with the local (7-oscillator 
algebra assigned to each site i = 1,. . . ,L. We take the 
quantum Lax operator in the form 



^^(«) = ( b, 1 



(2) 



which is equivalent to the Lax operator in [8,9,12]. The 
monodromy and transfer matrices are correspondingly 



T{u) = Ci{u)C2{u) ■ ■ ■ Cl{u) 



and t{u) = Trace T{u). The power-series expansion of 
the transfer matrix gives the commutative set of integrals 
of motion (IM), 



t{u) = 1 + (1 - q^)ur+ + u^-p|^^ + u^r^^^ 



(4) 



,2\„,L-l-i 



One additional integral of motion is given by the quantum 
determinant of the monodromy matrix: 



g-detT(u) = u^q^^ , M = . 

1=1 

The eigenvalue of M is the particle number N . 
The first non-trivial operators in (4) are 

1 ^ 
1 ^ 



(5) 



(6) 



where here and in the following periodic boundary condi- 
tions apply: £ = L + t. The Hamiltonian of the g-Bose gas 
is defined by 



(7) 



where /i is the chemical potential. Furthermore, the cur- 
rent operator J is given by J = ^ (P+ — V- ) and is 
equally conserved. 

^ Our definition corresponds to [4] by g — > b — » 

v/i-g-^B. 



The model can be solved by remarkably simple coor- 
dinate and algebraical Bethe Ansatze [8,9,12]. The most 
convenient form of the BA equations is the Baxter equa- 
tion relating an eigenvalue of the transfer-matrix and a 
polynomial Qiv of the power N , where N is the number 
of particles, 

t{u)QN{u) = QN{q^u) + u^q^''QN{q-^u) . (8) 

A simple derivation of the Baxter equation is given in the 
appendix. The roots of Qn are the Bethe roots. 



N 



Qn{u) = JJ(u-a;a) 



(9) 



0=1 



and the consistency condition for (8) gives the Bethe 
Ansatz equations 



quJa - q i^b 



q ^LOa - qOJb 



, n=l,...,N (10) 



(3) Set Wo = exp[zfca], then it follows from (7,8) 



N 



n = ^s{ka) , s{k) = - cos(A;) . (11) 



0=1 



2.2 Higher Integrals of motion. 

The power-series ansatz (4) has also been performed in [9] 
to obtain the Hamiltonian. To calculate all the other IM, 
it is convenient to rewrite the ^-operator as 

Ci{u) = uen + ub^iei2 + hee2i + 622, (12) 

where ea/3 is a 2 x 2-matrix with only non- vanishing entry 
1 in row a and column P, such that Cape-ys = SpjCoiS- 
Then one finds the coefficient of (w^~^) denoted as 

p|2) (-^(2)-) 



Vf^ = J2 (b£bt,+2 + b^b+^+ib,+2bt,+3) 
i=l 

p(2)= [p(2)jt_ 

(3) (3) 

Furthermore, the operators {T^- ) with coefficients 
{u^~^) read 

L 

■P+ ' = ^ (b^b^^+3 -I- b^b^^+ib£+2b^^+4 

e=i 

+ b£b^£+2b£+3b^i!+4 
+b£b^£+ib£+2b^£+3b£+4b^£+5) 

p(3) ^ ["p(3)"l 
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The operators v'^'' 



V. 



for general u are sums 



m2{e.) 



1=1 



(13) 



where the inner sum carries over all possible choices of 
pairs hn^i^e.)^'^ nk(t) such that v = md{i) — £+1 — Np, with 
Np being the number of bb^ -pairs. 

From Eq. (8), the eigenvalues of the operators are 
readily derived, namely for v = 2: 



and for v = 3: 

(1 - q^) 



2E 



E 



In the J^th order, one has all possible sums 



WoWo'Wa" 



(14) 



such that 1'j = y- 



3 Limiting cases and effective field theory 

In the limiting cases g ^ 1, g ^ 0, the Hamiltonian is 
written in terms of canonical Bosc operators, as shown 
in the first part of this section. In the second part, an 
effective field theory valid at low energies is given. 



3.1 Limiting cases 

To gain a first physical insight, we write the Hamiltonian 

in terms of bosonic operators a^, a]^ with [a^j On] = ^m,n- 
These operators act on states of J- as 

a\n) = \/n\n — 1), (^\n) = \/n-\- l|n + 1), 

a\n) ~ n\n). (15) 

On the other hand. 



n|n) = n|n), b|n) = \/\ — g^"|n — 1), 

bV) = + 

We can thus identify n = c^a. Let us set q = and 
expand 



(16) 



Then b, is written in terms of a, a}: 

OO 'OO 

b = ay^ X &„(2r/n)", b^ = ^ ^ 6„(2?yn)"a^ 
such that 

oc 

b^bt^+1 = 27] Y,^ebnbm{2vr+"' {neT (n^+i)'" aj+i. 

Following this scheme, we write down the Hamiltonian (7) 
explicitly for small r], i.e. q near one, including the order 



n = 



1 , 2 i 

2 — (^a^a^+i + a^a^-i J 

L 



(at 



0^ 



(17) 



7? 



2 r 
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+- (n^ + n|+i) aj+id^ + n£aj_,_ia£n£+i + /i.e. 



The first term in Eq. (17) is just the free hopping. The 
term linear in 77 describes assisted hopping processes: The 
amplitude is enhanced by the occupation number of the 
place to (or from) the hopping takes place. Note that 
this term acts repulsively. Such a term again is present 
in O {jf)- Additionally, hopping terms assisted with the 
square of the occupation number of the involved lattice 
sites occur. The sign of the 77^-term is negative, it lowers 
the energy. Generally, the interaction terms in order rf 
involve 2v + 2-many operators, which can be grouped into 
aj, a^+i and j/-many number operators n^, n^+i. Thus 
the interaction is local in each order and consists of these 
kinds of assisted nearest neighbour hoppings. 

Let us now consider the case where g > 0. Instead of 
(16) one expands 



Cnq 



n=0 



such that 



.n=0 



The first two orders in this g-expansion of the Hamiltonian 
read 



1 



1 



alae+i 



\/ni+T 



n=0 



(18) 
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The next term stemming from the "Pi-operators is of the 
form 



1 ^ 



1 



H —a^ae+i 



1 



The operator (18) ck-scribcs a bosonic model with nearest- 
neighbour hopping where the hopping amphtude is in- 
versely proportional to the square-roots of the occupancy 
of the lattice sites involved. This means that the hopping 
amplitude equals one regardless of the occupancy. 



3.2 The effective field theory 

It is instructive to consider the Hamiltonian (7) within the 
density phase representation of the operators b, b^. We 
define operators E^^^ by [10,20] 

OO GO 

£;(+) :=^|n+l)(n|, E^-'> :=^|n)(n + l|, 



n=0 



n=0 



such that at = E'^+^^/n+l, a = ^/rT+TE^f-) and 



E^^\ However, the E^"^^ are not 
= 1 , = 1 - |0)(0|. (19) 



Note that (E(±))t 
unitary: 



Then 



b = Vl-g2(n+l)£;(-), bt = E(+)V1-9^("+^^- 

Let us now formulate an effective field theory at small 
energies, i.e. at long wavelengths and small rj. Therefore 
we pass to the continuum via 



X = I ■ A, ni ^ A- p{x), 
[p{x),E±{y)]^±5{x-y)E''^\y). 



(20) 



In order to be able to treat E± as unitary operators, we 
set 

p{x) ~ n + 4>{x) 

with {4>{x)) <C n and {dx4>{x)) <C 1, such that (p en- 
codes the small-amplitude, long-wavelength-oscillations of 

the density around its equilibrium value n. Here n should 
not be confused with the label of the Fock states. Thus 
we neglect the vacuum component in (19) and set E± — 
g=Fi^r(a:)_ Because of the commutation relation (20), we 
have [p{x),n{y)] = iS{x — y). This leads us to 

b = (l - ^) e'" + O {a^'^) 

bt = e-'" (l - ^) + O {A''^ 



We now expand —II{x) + II{x + A) = AdxU and discard 
higher-order terms, since in the continuum limit A 
they do not contribute. We end up with 



n = - 



A , 
-n(5.iT)' 



-r]{n'^ + 2n(t){x) + (t)^{x)))dx. (21) 

The first and the third term yield the ground state energy 
density e = —n + ijn?. It will be discarded in the rest 
of this section. Since 4'{x) is oscillatory, J (j){x)dx = 0. 
Setting 7] = cA, we end up with 



n 



A 
"2 



/ 



n{dxn{x)f + 2c(j)'^{x) 



dx. 



(22) 



Before proceeding further, let us note an alternative way of 

carrying out the continuum limit of the Hamiltonian (17): 
Set ae = VA^{x), where a; = £ • Z\ as in Eq. (20). Then 
the field operators ^, obey the commutation relation 

[^(x),W^y)] = 5{x-y). If we expand a^+i/v^ = 'P{x) + 

Adx'l'{x) + ^dl\l/{x) + . . ., and again put rj = cA, then 
after an integration by parts and within the same order 
as (22), 

n = ^ I [dx^\x)dxnx) 



The corresponding quantum 
Hamiltonian Hn reads [16] 



+2c'^\x)^\x)\^{x)'^{x)\ dx. (23) 
mechanical AT-particle 



n 



N 



i=l l<j<k<N 



N 



(24) 



In [14,11] the equivalence between (22) and (23) has been 
stated. Thus the Bose-gas in the continuum with small 
coupling constant presents an effective low-energy model 
for the Q'-deformed Bose-lattice model. 

We continue the investigation of (22). For the case of 
notation, we set A — I and concentrate on the Hamil- 
tonian density H, defined by H =: / if da;. To bring the 
model (22) into the standard form, we define the charge 
velocity 



and the Luttinger parameter 
K = 



(25) 



(26) 



Then, by scaling 77 = II' /Vk, (j) = \/K<j)' , we obtain 77 = 
^^ (^{d^n' (x))^ + {(j)' {x))^'^ , which is very closely related 
to the Gaussian model. We will comment on it further 
below. Including a chemical potential term — mX]^=i'^^ 
adds a contribution — /u(n -|- \fK4>') to 77. Upon shifting 
6' ^ 6' + , we arrive at 



H=^(^{dxn'{x)f + {cP'{x)f) 



2vc ' 



(27) 
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such that for the charge susceptibiHty one obtains Xc = ^• 
In the low-couphng hmit considered here, 



Xc 



1 

2r]' 



(28) 



Also note that from (27) it follows that the specific heat at 
fow temperatures behaves as limT->o C{T)/T = Let 
us emphasize again that the results (25,26,28) arc valid in 
the limit ^ 0. Higher-order contributions can be princi- 
pally obtained by pursuing the field-theoretical approach 
to next-leading operators. However, it is much more con- 
venient to make use of the integrabihty of the Hamiltonian 
(7). This will be done in the next section. 

In the last part of this section, the low-temperature- 
asymptotics of the (/-deformed Boson-model will be con- 
sidered in the framework of the effective low-energy model. 
Therefore it is helpful to introduce chiral components 
4>R,L- We shortly remember how this is done for the Gaus- 
sian model, defined by 

[<i>G{x), ndy)] = iS{x-y). 

Here it is convenient to introduce a field Oq such that 
n^G = dxOc- To diagonalize the Hamiltonian and to obey 
the commutator, one performs a mode expansion of the 
fields (pG-.^Gi [19]. As a result, 9g = (j^L — 4>Rj 4>G = 't'L + 
(pu, where <pR {4>l) are the chiral components of the Bose 
field (f)G (i.e., the right- and left-moving parts). 

In a very similar manner, the model (27) can be 
treated. Namely, the diagonalization and the commuta- 
tor [(j)' [x) , n' {y)] = iS{x — y) lead to a mode expansion 
which now suggests to introduce 9' such that </>' = dx0' 
and 6' = (f)L + (t>R, n' = (f)L — (f>R- Prom this it follows that 

= !f = e'(^--'^«)/v^(n + V1^(0L + 0ij))i/2 (29) 

^=^^ = {n + VK{cPl + </)fl)) V^e-'^-^--^^)/^. (30) 
V A 



Based on these operators, the effective field theory in the 
continuum (that is, for low energies) is formulated. In the 
continuum limit, all operators that are not forbidden by 
symmetry are allowed - the model (27) only constitutes the 
leading order. To gain insight into the low-temperature be- 
haviour of the specific heat beyond the leading order, one 
has to consider next-leading contributions to (27). This 
has been done for the Heisenberg chain in [19,21]. Since 
the low-energy effective theory of the Heisenberg chain is 
a Gaussian model with next-leading terms as well, similar 
considerations apply in the present case. Let us enumer- 
ate the possible allowed operators. In the leading order 
(i.e. with scaling dimension 2) these are (dx4'L,R.)'^ ■ Com- 
binations of uneven numbers of dx<t>R.L a-i'c not allowed, 
so there are next-leading operators of scaling dimension 

4: [dxc^LM)'.. {dxH)yM\ {dl4>L.R)\ {d^MKd'M- 

Other relativistically invariant operators are cos a((j)ii + 
4>i,) and cosa{4>R — 4>l)- The first ones cannot appear be- 
cause it is impossible to built them from (29), (30). The 



latter ones are forbidden by the J7(l)-symmetry of the 
model: namely, multiplying a {a^) by e'''' (e"'''') does not 
change the model. Prom (29), (30) one sees that such a 
multiplication is equivalent to shifting (p^^ (j)R in opposite 
directions. However, an operator cos «(</>/{ — 4>l) does not 
preserve this symmetry. We conclude that the next-leading 
operators are of scaling dimension 2n, n > 2. According 
to the calculations done for the Heisenberg chain in [19, 
21], the low-temperature expansions of the specific heat 
and charge susceptibility then read 



C(T) 



K 



X. 



.(r) = ^5,r2^Jo = -. 



The constant can be calculated from the low- 
temperature expansion of the thermodynamical poten- 
tial g at constant n, g{T) + jjn 
X-^{T) = dlg{T), one obtains 



e - ^T^. Since 



Si 



K 



(31) 



The Bethe Ansatz solution allows it to verify numerically 

that 

C{T)/T = + 7ir2 ^ ^^(y^ ^^^2 ^ ^32) 

3Vc Vc 

with unknown coefficients 71, 62- 

4 Ground-state properties and 
thermodynamics from Bethe ansatz 

This section contains two parts. In the first part, ground- 
state properties are considered, whereas the second part 
focuses on the thermodynamics. 

4.1 Ground-state properties 

To investigate the ground-state properties, we start from 
the BA-equations (10). We rewrite them using the param- 
eterization LUa = exp [ika], q = exp [—r]]: 



JkaL 



n 



3in(^ 



\ sin (^ii_A _ i^) 



a=l,...,7V. (33) 



Using the same argument as in [16] (chapter 1.2). one can 
prove that ka € M, a = 1, . . . ,N. We follow the common 
procedure of introducing the density p{k) of roots in the 
thermodynamic limit 

Si I I I [ I A 7 / 

-p{q)dq (34) 



sinh 2i] 



2tt 2'k J_g cosh2r] — cos{k — q) 



N 

-=:n 

E 



r p{k) 

J-B 



dk 



cos k p{k) dk . 



(35) 

(36) 
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Let us consider the limiting cases q ^ and q 
q ^ 0, that is rj — > oo, 

= i!^ + 0(^2^ cos fc), 



1. For dressed energy, parametrized by the quasimomentum k, is 
£o(^) = ~ cos k — 11 + 

sinh 2r} 



n + 1 
n+ 1 



O (q'n) 



sm + O iq^n) 

TT n+1 ^ ^ 



(37) 
(38) 



For g = 0, the root density is constant, a property known 
from the free fermion gas, and the energy is bounded with 
respect to the density. The order ~ q^n is consistent with 
the corresponding Hamiltonian (18). Note that for q = 
0, the Bethe ansatz equations (33) are solved explicitly. 
Namely, for g = 0, these are solved by (for the ground 
state) 



N + L 



{-N + 2a-l), a = l,...,iV. 



This result can also be interpreted as if the particles were 
free (i.e. fc^ — 2'na/ L), and each mode has occupancy num- 
ber n+l. The above q = 0-results have already been found 
in [5,9]. There it was also shown that one can solve the 
BA equations for g = not only for the ground state, but 
for general excited states [7, 8] . This allows it to calculate 
correlation functions exactly [7,8]. 

The opposite limit g — > 1, that is ?7 ^ 0, is technically 
more involved. The integration kernel in Eq. (34) becomes 
singular in this limit. Instead of dealing with the integral 
equation (34), it is more convenient to perform a small- 
expansion of Eqs. (33). This technique has been applied 
successfully to quantum gases in [13,1-3]. Here we find 
that 



(39) 



Thus 



where the qaS satisfy the equation qa = X^^^o ga-gh ' 
the qaS are the roots of the Hcrmite polynomial Hj\f{q). 
Inserting Eq. (39) into Eq. (35) yields 



e = —n + r]n 



(40) 



which agrees with the field-theoretical result for weak cou- 
pling, stated after Eq. (21). It also coincides with the 
mean- field approximation of the Hamiltonian (17) includ- 
ing the order linear in rj. Furthermore, from the proper- 
ties of the gaS, the root density behaves semi-circular-like: 
p{k) 



\J 8r]n — fc^ . Note that essentially the same be- 
haviour has been found for the repulsive Bose gas [13], 
again implying the continuum limit (24) for g ^ 1. 

The model (7) is conformally invariant at < g < 1 
(for q = 1, the dispersion relation is quadratic and thus 
conformal invariance is broken). This can be seen most 
easily from the BA solution: The low-lying elementary ex- 
citations arc the addition or removal of a particle with 
momcintum around the two Fermi-points ±Pf with pp = 
with a linear dispersion relation. The corresponding 



7rn 
2 



27r/ 



27r J_g cosh2?7 — cos(fc — q) 



Soiq)dq, (41) 



where /i is to be determined such that eo(S) = 0. The 
velocity Vc of these charge excitations (which is identical 
to the sound velocity) reads 



Vc 



" 



dp 



P=PF 



dkp{k) 



k=B 



2-Kp{k) 



k=B 



After the first equality sign, the usual relation to obtain 
the velocity for massless excitations with momentum p 
and energy Sq has been used. Then p and Sq have been 
parametrized by the quasimomentum k, and finally the 
relation between p and k has been inserted. In analogy 
to the consideration of limiting cases in the ground-state 
energy, we obtain for Vc- 



1 



1 



fc(n)|g- 



The last equation has been obtained from the equation for 
the sound velocity [17] 



Vc = 



mn 



1/2 



(43) 



where the mass m = 1 here and the pressure P — —OlE, 
with E — eL the total ground state energy. Thus inserting 
Eq. (40) into Eq. (43), the second result in (42) is obtained. 
Note that it agrees with Eq. (25), derived in the field- 
theoretical approach. 

The charge velocity is inversely proportional to the 
charge susceptibility Xc at T = 0. To see this, we cite the 
result originally derived for the one- dimensional Bose gas 
[16], but also vahd here 



Xc 



irvc 



(44) 



where Z = ^{B) and the function ^{k) = 2iTp{k). The 
Luttinger parameter K is thus obtained as 

K = ZVtt = VcXc- (45) 

In the special case g = one has Z = n + 1, so that 



Xc(g = 0) 



TT sm 



(46) 



which diverges ~ l/(n7r2) at n ^ and as ~ ti^/tt^ at 
n oo. 

To obtain a qualitative picture for the various quanti- 
ties at arbitrary < g < 1, Eq. (34), and the derivative of 
Eq. (41) with respect to the quasimomentum are solved 
numerically. The results are shown in Figs. 1, 2, together 
with the field-theoretical results in the weak-coupling limit 
g^l. 
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Fig. 1. a) The ground state energy e{q) for different densities. 
The straight lines are the q 1 and g — > approximations, 
respectively (cf. Eqs. (38), (40); b) the charge susceptibility 
Xc{q), Elq. (44). The thick blue line near q = 1 shows the field- 
theoretical result (28) which is independent of n. The straight 
lines near q = depict Eq. (46) . 





Fig. 2. a) The charge velocity Vc{q) for different densities 
n. Additionally, near q = 1,0, the analytical results (42) are 
shown; b) the Luttinger parameter Kc- The lines near q — 1 
denote the result (26) from field theory, whereas near q = 0, 



the values for K\ 



q=0 



(n + 1) /tt are shown. 



We finally comment on the parameters n, fi, B. Since 
dBTi = 2p{B)Z > 0, the particle density is a mono- 
tonically increasing function of B. As p{k),e{k) are 2tt- 
periodic, the integration boundary B is restricted to i? G 
[0,7r]. If B = 0, one has n = 0. For B = tt, Eq. (34) 
together with (35) is satisfied only for n — s- oo. Let us 
now determine the corresponding /i-values. Thermody- 
namic stability requires S^n > 0. As stated under Eq. (41), 
the condition e{B) — fixes fi. Thus ^ = —1 for B = 0. 
If B = TT, Eq. (41) can be treated by Fourier-transform, 
yielding fi — 0. These findings are summarized in table 
1. In the special case q = 0, we have from Eq. (37) that 



n 





oo 




-1 





B 





TT 



Table 1. The ranges of the parameters n, fj,, B (at T = 0). 



B = rnr/ {n+ 1), and from Eq. (38) one derives p, = — 9„e, 



li(n) = — [ — sm 



1 mr 
■ cos ■ 



TT n + 1 n + 1 n + 1 



4.2 Thermodynamics 

As stated in the previous section, all solutions to Eqs. (33) 
are real. This means that in order to describe equilibrium 
thermodynamics in the thermodynamic limit, we can fol- 
low the Yang- Yang approach [22] which has been devel- 
oped to calculate the thermodynamical potential g{T,ii) 
for the one-dimensional Bose gas (24) . Here T is the tem- 
perature and ^ the chemical potential. As a result, one 
obtains 



l{T,p) = - 



t£,„(i+c-*V.) |, 



where the function e{k) is obtained from 
e(fc) = — cos k — ji 

-T 



(47) 



sinh 277 



, -lnfl + e-(^)/^) ^, 
^ cosh 2?7 — cos(k — q) V / 27r 



In the zero-temperature limit T — > 0, Eq. (47) turns into 
Eq. (41), with limr^o £ = £o- 

We are interested in thermodynamic quantities like 
the specific heat C(T) and the charge susceptibility Xc[T) 
(both per lattice site) at fixed particle density n. Let us 
first consider C{T). Within the grand canonical ensemble 
considered here, it is calculated using thermodynamic re- 
lations, cf., for example, [15]. From the thermodynamical 
potential g(T,ij), one obtains the entropy per site S at 
fixed /i and the particle density n 



S = drgl 



To keep n fixed, the chemical potential acquires a T- 
dependence, p = p{T), such that 



drfJ-ln 



df^n\T 



In the numerics, the partial derivatives are calculated by 
setting up the corresponding integral equations and iter- 
ating these, using the Fast-Fourier- Transform (FFT) al- 
gorithm. 
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We have performed several consistency checks of the 
results following from the iteration procedure in certain 
limiting cases. First of all, at low temperatures the nu- 
merical data are compared with the analytical result 

C(T) = ^T, (48) 

which is derived from a saddle-point integration in (47) at 
T ^ 0. Note that this result agrees with the prediction 
from conformal field theory, given that the system is con- 
formally invariant at a finite Vc (that is, at < 9 < 1). 
Furthermore, at q — 0, 

J-7T ^ ^ 27r 

such that the thermodynamical potential is given by 

g{T,n)^-^J ln(l + e^™"'=+'5'^-'^s)dfc. 

These equations are numerically easier to solve because 
no convolution has to be done. Especially, in the high- 
temperature limit T — !■ cx) and large densities n 3> 1 

9-.r(-i„„-i + o(-l))-f^ (49) 

Remarkably, the entropy ^ Inn and the specific heat 
C ^ 1/(4T^) in this limit. From the numerics for gen- 
eral q > 0, we find that C ~ a/T'^ with an interaction- 
dependent constant a. 

As an illustration. Fig. 3 shows the specific heat at 
the values g = and q = exp [—0.5] for different densi- 
ties. It is interesting to note that at high densities, the 
specific heat displays a double-peak structure. We inter- 
pret this behaviour as follows. At lowest temperatures, 
the particle-hole excitations around the two Fermi points 
are activated, resulting in the slope oc T of the specific 
heat. At higher temperatures, all particles at each lat- 
tice site are thermally activated. The activation of this 
"bulk" causes the second peak, which is clearly discernible 
at high densities. The numerical accuracy of our data is 
sufficient to confirm Eq. (32). The next-leading term oc 
in C{T) - Ttt/{3vc) is shown in Fig. 4. 

We now turn our attention to the charge susceptibility 
Xc = ~d^g{T, Ij), calculated at fixed n. Again for general 
q, the corresponding integral equations are solved numeri- 
cally by iteration, using the FFT. The T = 0-limit is given 
in Eq. (44). For g = at high temperatures, one obtains 
from Eq. (49) that Xc ^ r? jT . In Fig. (5), the charge sus- 
ceptibility is shown for q — Q and q = exp [—0.5] at differ- 
ent densities. The leading T-dependent contributions to Xc 
are depicted in Fig. 6. They confirm the field-theoretical 
prediction Xc(r) — Xc(0) = b\T'^ at low T. Especially, for 
g = 0, the constant given in Eq. (31), is calculated 
explicitly from Eqs. (42), (46). This allows it to confirm 
numerically Xc{T) - Xc(0) = (JiT^ + b^T^ for 9 = 0, with 
an unknown coefficient 




T T 



Fig. 3. The specific heat C(T) at a) g = 
and b) g = exp [—0.5] for difi'erent densities 
(in a): n = 0.00774, 0.649, 2.09, 2.95, 5.86, 13.9 
from bottom to top near T = 0; in b): n = 
0.01479, 0.2306, 1.059, 2.356, 3.441, 4.658, 8.399 from bottom 
to top near T = 0. The slope at T > agrees with (48). 




T T 



Fig. 4. Next-leading T-contribution to specific heat: C(T) — 
r (X a) g = and b) g = exp [—0.5]. For comparison, 

the dashed red lines are 1000 ■ in a) (100 • in b)). This 
confirms that the exponent is 3 in all cases. 



5 Conclusion 

We investigated the exact solution of the q-deformed Bose 
gas. Emphasis was put on the higher integrals of motion 
by using the algebraic BA. We furthermore presented a 
detailed derivation of the effective field theory for low en- 
ergies, with comparisons between the field-theoretical pre- 
dictions and the BA. These comparisons were done both 
for ground-state properties and for thermodynamic quan- 
tities. Open questions include, for example, the calculation 
of the coefficients of the leading irrelevant operators, a 
high-temperature expansion of the thermodynamic quan- 
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0.5 

T 



Fig. 5. The susceptibility Xc at a) g = and h) q — exp [—0.5] 
for densities n — 0.5, 1, 2, 4, 7 from bottom to top. Diamonds 
on the t/-axis denote the T = 0- values obtained from (44). 
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Fig. 6. Leading T-contribution to the susceptibility Xc(T) — 
Xc(0) « at a) g = and \>) q — exp [—0.5]. In a), the 
blue dotted lines are the exact result Eq. (31), obtained from 
Eqs. (42), (46). The inset shows Xc{T) - Xc(0) - ^iT^. For 
comparison, the dashed red line is lO** • showing that the 
exponent is 4 in all cases. The red dashed line in b) is 10 ■ T^, 
confirming the exponent 2. 



titles and the problem of formulating a quantum transfer 
matrix for a lattice boson model. 
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A Derivation of Baxter equation 

Consider the following systems of equations for vectors 

!</.), |<^'), 10") e^: 



u uW 
b 1 



i') (/?', 1) , 



b 1 



1 

-13' 



Solutions of these systems are 

oo 

10) = \q^p,p>{u)) ^/3'"c„ 



1 

-/3 



13 

r 



(50) 



bt"|0) (51) 



and 

10') = \^f3,P'iq'u)) , 10") = V"l0/3,/3'(9"'«)> , (52) 

where c„(u) is defined by the recurrent relation 

(l-q^")c„(M) = [1 - u)Cn-l{u) + UCn-2{u) (53) 

with initial conditions cq = 1, c_i = 0. Note that c„(u) is 
an nth power polynomial of u. 
Let further 



l^piu)) = |0/3i,/32(w)®'^/32,/33 (")®- 

Since 

1 ^ (?)(/3i,l) + 

then 

Trace r(u) = (/3i, l)r(u) 
and it follows from (50) and (54) 



(54) 



1 

-A 



(1,0) 



(55) 



1 



(56) 



t{u)\-Pp{u)) = \<Pp{q'u))+u\^-^\<Pp{q--'u)) . (57) 

Let further (!f'Ar| be an eigenvector of t{u), N being the 
number of bosons. Set 



Qi3m{u) = (!fjv|<?/3(u)) 



(58) 



Since the set of (3i, P2, ■ ■ ■ , Pl is generic, Qi3m{u) 7^ 0. 
The coefficients Cn{u) are polynomials, therefore Qfs^Niu) 
is a well defined polynomial of u of the power N. It sat- 
isfies the Baxter equation (8). Obviously, the set of /3 in 
Qi3.n{u) corresponds to a spectral parameter independent 
normalization factor: 



Qp,n{u) = KpQn{u) 



(59) 
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